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Observation of strong final-state effects in tt + production 
in pp collisions at 400 MeV 
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Differential cross sections of the reactions pp — ► dir + and pp — ► pn7r + have been 
measured at T p = 400 MeV by detecting the charged ejectiles in the angular range 
4° < &Lab < 21°. The deduced total cross sections agree well with those published 
previously for neighbouring energies. The invariant mass spectra are observed to 
be strongly affected by A production and NN final-state interaction. The data are 
well described by Monte Carlo simulations including both these effects. The ratio 
of pp — > pnir + and pp — > dir + cross sections also compares favourably to a recent 
theoretical prediction which suggests a dominance of np-production in the relative 
3 S'i-state. 
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The single pion production has received renewed interest in recent years both from ex- 
perimental and theoretical points of view. Theoretical aspects currently under discussion 
include possible heavy meson exchange, the nature of the ttNN vertex and the role of final- 
state interactions (FSI) in these reactions. Experimentally the availability of storage rings 
has opened the possibility to measure single pion production with high statistics even close 
to threshold fl|j2|,[l,iJ • These data show an energy dependence of the total cross section near 
threshold, which deviates substantially from phase space suggesting a strong influence of the 
NN FSI. In this Letter we show that FSI effects can be explicitely observed and identified in 
the invariant (M) and missing mass (MM) spectra of the reaction pp —>■ pnir + . The energy 
of T p = 400 MeV is already high enough for A-production to be observed clearly in M nn + 
and M p7r +, whereas the pn FSI still strongly influences M pn even though this energy is well 
above threshold. 

The measurements have been performed at the CELSIUS storage ring at T p = 400 MeV 
using the WASA/PROMICE detector setup including a hydrogen cluster jet target. Details 
of the detector and its performance are given in ||. For the data presented here only the 
forward detector has been utilized, which allows the determination of the four-momentum of 
charged particles in the angular range of 4° < ®Lab < 21°. It is composed of a tracker with 
proportional counter straw chambers for an accurate determination of particle trajectories, 
followed by segmented scintillator trigger and range hodoscopes for dE and E determinations, 
respectively. Particle identification has been made by use of the dE-E method. Fig. [I] shows 
a three-dimensional plot of the dE-E spectrum. Deuterons, protons and pions appear well 
separated. For the 7r + identification the delayed pulse technique (observation of the delayed 
pulse from fi + decay following the 7r + decay) has been utilized in addition. Though protons 
and deuterons appear to be well separated in Fig. [I], their separation is not perfect and 
the contamination of reconstructed pnir + events with d,7r + events, and vice versa, is non- 
negligible. In order to get rid of this contamination we have requested the dn + events to 
be planar (A$ = 170° — 190°) and the p7i + events to be non-planar. This way mutual 
contaminations could be kept below 1%. The final number of good events which passed all 
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criteria has been about 10 5 for dir + and 7 * 10 4 for pn7c + . The integral luminosity has been 
determined to better than 5% by the simultaneous measurement of pp elastic scattering 
and its comparison to literature values 0. Detector response and efficiencies have been 
determined by Monte-Carlo (MC) simulations, utilizing the program package GEANT || 
including the treatment of secondary interactions in the detector. 

The data obtained for pp — > dn + are shown in Fig. ^|. In the upper part the 7r + missing 
mass MM n + which gives a peak at the position of the deuteron mass is displayed together 
with the corresponding MC simulations. The good agreement demonstrates that the detector 
response is understood to high precision. The measured tt + angular distribution in the 
center-of-mass (cm.) system is displayed in the lower part of Fig. ^together with the SAID 
phase shift prediction [0. Note that the angular distribution is symmetric about 90° due to 
the symmetry in the entrance channel. The limited angular range of the data results from 
the experimental requirements of 4° < Qiab < 21° for both d and tt + . The experimental 
points are compatible with SAID, and we may use the angular dependence of the latter 
to extrapolate our data to Att. This way we obtain a total cross section a(pp — ► d% + ) = 
0.78 mb, which compares quite favorably with the SAID value of 0.82 mb. Whereas the 
statistical uncertainty is less than 1%, the systematic uncertainty is estimated to be about 
7% comprising uncertainties both from the determination of the luminosity and from the 
handling of deuteron breakup in the detector by the MC simulation. 

Results for the pp — > pnn + reaction are comprised in Figs. |3|-|5] . The invariant and 
missing mass spectra reconstructed from the measured pir + events are shown in Fig. |3| 
together with curves from MC simulations assuming pure phase space (dotted) and including 
either A ++ (A + ) excitation in the p7r + (nn + ) system (dash-dotted) or the s-wave FSI in the 
pn system (dashed). For the latter the effective range approximation of Migdal- Watson 
type as given in || has been used. The only free parameter, R, accounts for the size of the 
interaction region, from which the two nucleons emerge, and effectively subsumes also the 
part of the reaction, where the two nucleons are not in relative s-waves. For a point-like 
vertex one finds R ~ 0.8 fm ||. Here in the MC simulations R has been adjusted for best 
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reproduction of the data resulting in R = 2.5 fm, if we assume the s-wave pn system to be in 
pure 3 Si states as previous work suggests [@,|nj. This assumption will also be corroborated 
by a comparison of the dir + and pnir + channels below. The value for R is compatible with 
those obtained in analyses of the NN FSI in the deuteron breakup on the proton at low 



energies ||11|| . The A excitation in the exit channel is assumed to take place between the 
7T + and either the proton or the neutron according to the isospin ratio of 9 : 1 with the ir + 
being in relative p-state in the ttN system. Note that parity conservation requires the pion 
also to be in a p-wave relative to 3 Si pn final state. 

The MC simulation including both FSI effects together is displayed in Fig. |3| by the 
shaded histograms which reproduce the data very well. The invariant mass spectrum M p7r + 
and the missing mass spectrum MM P (corresponding to M n7T +) are peaked towards large 
masses, whereas MM 7r +( corresponding to M pn ) is strongly peaked towards low masses in 
contrast to phase space distributions (dotted lines). In all three spectra both FSI effects play 
a significant role. However, the fact that the peaking towards large masses is much more 
pronounced in M pn + than in MM p , clearly exhibits a strong influence of the A production, 
since it is nine times more likely in the pn + system than in the mr + system. This observation 
in M p7T + and M n7T + is strongly different from the one at T p = 310 MeV, where both spectra 
have been observed to be of comparable shape However, it is in agreement with the trend 
observed for T p < 330 MeV || . From the partial- wave analysis of those data it was found 
that the resonant p-wave contribution, though still small at these energies, is steeply rising 
with incident energy. From this analysis the A contribution is expected to get dominant 
at energies of T p pa 400 MeV (see Fig. 20 of Ref. [fj]). This is indeed what we observe in 
our data. We note that the strong influence of the FSI on the pion energy distribution is 
already apparent in the dE-E plot in Fig. [I], where the pions of the pnn + reaction are seen 
to be concentrated towards the peak of the d,7r + reaction. 

In Fig. ^ the experimental pion and proton angular distributions, corrected for detector 
efficiency and acceptance, are compared to MC simulations for the reaction pp — > pmi + . 
As in pp — > dir + the angular distributions have to be symmetric about 90°. The proton 
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angular distribution is essentially isotropic being affected only slightly by the A excitation. 
The pion angular distribution, on the other hand, depends strongly on the A excitation. 
Conventionally the pion angular distributions in single pion production are parametrized by 
a(0 7r ) ~ 1/3 + 6cos 2 7r , where b is the so-called anisotropy parameter. Previous analyses 
yielded values for b up to 0.5 for np — > nn-K + ]I2] and 0.3 for pp — > ppn° |L3| with the 
maximum b being reached near T p m 550 MeV. The b values are observed to decrease with 
decreasing T p . At T p = 460 MeV, the lowest energy analyzed in those studies, b gets as small 
as 0.1. On the other hand the IUCF measurements [p],[| of pp — > pmr + show a strong rise of 
b already close to threshold reaching b = 0.23(6) at their highest energy of T p = 330 MeV. 
Whereas from the latter we would expect to find already a quite substantial value of b at 
400 MeV, the previous analyses would suggest rather a small value. Hence we show in Fig. |] 
two MC calculations including pn FSI and A excitation, one with b = 0.1 (dashed lines) 
and one with b = 0.4 (solid) — in addition to the phase-space expectation (dotted line). In 
the measured range of 0^ the data clearly prefer the larger b value. From the cos 2 (0„-) plot 
of the pion angular distribution in Fig. |] it is readily seen that b = 0.7 would be an upper 
limit compatible with our data. For a more precise determination of b larger pion angles 
are necessary, which are not covered by this measurement. Hence also the determination of 
the total cross section from our data is not independent of the assumption for b. Whereas 
for b = 0.1 one would obtain a value of a(pp — > pmr + ) = 0.73(4) mb, our data favour values 
of 0.62(4) and 0.57(4) mb for b = 0.4 and 0.7, respectively, where the assigned uncertainty 
is essentially due to that of the luminosity. The latter values are in good agreement with 
literature values at neighbouring energies (see |L|). 

The close relation between o?7r + channel and pn7r + channel as the breakup channel of the 
former has recently been pointed out by Boudard, Faldt and Wilkin ||. Assuming that the 
final pn system is in the 3 Si state, and that the pion production operator is of short range, 
they derive a simple relation (eq. 6 of ref. ||) for the ratio of differential cross sections 
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defined by 



R(Q) = 2nB d 



dVt^dQ 



(pp —>■ prm + ) j 



da 



da 



7T 



(pp — > diT + ) 



p(x) sfx 
p(— 1) 1 + X 



which should be independent of the pion scattering angle. Here Bd denotes the deuteron 
binding energy and Q = M pn —md the excitation energy in the pn system, x = Q/ Bd and p(x) 
is the pion cm. momentum. Our experimental result for R(Q) is plotted in Fig. || together 
with the prediction of Ref. ||. Though there is good qualitative agreement, the data exhibit 
a significantly steeper slope than the prediction. In particular, the data are higher at low 
Q, where the approximations made in Ref. U should be valid best. Experimentally the 
determination of R(Q) from the simultaneous measurement of both reactions is expected 
to be particularly reliable, since uncertainties in the determination of the luminosity and 
detector response cancel to large degree in the ratio, the only major source of possible error 
being the treatment of the deuteron breakup in the detector. However, since our result for 
a(pp — > dn + ) agrees well with the literature values at neighbouring energies, we do not see 
a significant problem there either. A closer inspection of Fig. 1 in ref. |5) indicates that 
experimentally R(Q) is not fully angle independent. The TRIUMF data [[KJ plotted there 
exhibit a systematic trend around the maximum of R(Q): the data taken at small angles 
(@La& = 46°, 56°) lie systematically above the theoretical curve, whereas the ones taken at 
larger angles (73° — 88°) lie significantly below. Our result for the forward angular range fits 
very well into this trend in the TRIUMF data. Since in the calculation of R(Q) only the 
isoscalar 3 Si channel of the pn system is considered, the surplus of approximately 10% in 
our data at small Q may be associated with contributions from other partial wave channels, 
notably the 1 So channel of the pn system. This conclusion conforms quite nicely with the 
expectations from the partial- wave analysis at lower energies 0. The observed discrepancy 
between calculation and data at higher Q may be associated with effects from higher partial 
waves as well as with kinematic approximations made in the derivation of the theoretical 
expression for R(Q) ||. 

Summarizing we observe in the exclusive measurement of pp — > pmr + at T p = 400 MeV 
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strong FSI effects in the invariant and missing mass spectra which are identified as being 
due to pn FSI and ttNA excitation. The dominance of the latter is in agreement with 
expectations from partial-wave analyses close to threshold. The influence of the former is 
in agreement with expectations for the Migdal- Watson effect in the pn system, if s-waves 
are dominating there. The observed ratio R(Q) for the exit channels pnir + and dn + is 
compared with the prediction of ref. || relating the channel of the bound system with its 
breakup channel. Although the general agreement is good, there are significant differences 
which call for a more refined theoretical treatment of both channels. 

We gratefully acknowledge valuable discussions with Colin Wilkin. We are also grateful 
to the personnel at the The Svedberg Laboratory for their help during the course of this 
work. 
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FIGURES 

FIG. 1. Plot of the dE-E spectrum for two charged particles in the forward detector (FD). The 
energy loss dE has been measured by the forward trigger hodoscope (FTH), see ref. ||. Energies 
are given in units of GeV. 

FIG. 2. Top: Spectrum of the ir + missing mass MM W + as obtained from identified dir + events 
together with the corresponding MC simulation (shaded area) of the detector response. Bottom: 
Measured tt + angular distribution for pp — > dn + in comparison with the SAID phase shift prediction 

FIG. 3. Spectra for pir + invariant mass M p7T + and missing masses MM p and MM W + for 
pp —* pnir + reconstructed from the pir + events detected for 4° < @Lab < 21°. Corresponding 
MC simulations are shown assuming pure phase space (dotted), including either pn FSI (dashed) 
or A excitation (dash-dotted) and both together (shaded histograms). 

FIG. 4. Angular distributions of protons (top) and pions (bottom) for pp — > pnir + reconstructed 
from the detected p7r + events. They are plotted in dependence of cos 2 (0) and are compared to 
MC simulations assuming pure phase space (dotted) and including both pn FSI and A excitation 
with b = 0.1 (dashed) or b = 0.4 (solid). All simulations are normalized to an integral cross section 
of a = 0.62 mb. 

FIG. 5. Ratio R(Q) of the pp — > pnir + and pp — * dir + differential cross sections. The solid 
curve shows the prediction of Ref. ||] . 
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